Persistent current on a disordered mesoscopic ring with an embedded quantum dot 
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We calculate the ensemble averaged persistent current on disordered mesoscopic rings with an 
embedded quantum dot. We model the quantum dot as a single resonance and use Random Matrix 
Theory to model the impurities in the ring. Using Efetov's supersymmetry technique, we develop an 
analytical expression for the current. We find not only one, but two resonance peaks in the current. 
This is interpreted as quantum interference phenomenon. 



INTRODUCTION 

Given the rapid development of the technology in the 
fabrication and manipulation of semiconductor struc- 
tures, there has been great interest in the study of dif- 
ferent mesoscopic systems, specially involving quantum 
dots. Many atomic effects as the Kondo and Fano effects 
can also be observed in these systems and in controllable 
way lj. Those work deal mainly with the well known 
Aharonov-Bohm interferometer, which is basically a ring 
attached to at least two external leads. Then it is possi- 
ble to measure, for instance, the interference pattern of 
the electronic wave function as a function of an external 
magnetic field, or the dephasing of this very function by 
the tunneling of the electron through the quantum dot, 
among many other features. 

We are here basically interested in an isolated ring with 
an embedded quantum dot, and its effect on thepersis- 
tent current. Some recent works S 0, H B H @ 
have driven attention to this problem, specially in the 
Kondo regime. However none of them consider the dis- 
order in the ring, and are solved numerically by means 
of the slave-boson mean-field theory. We propose a new 
approach to the problem, and also find some interesting 
typically quantum features of the current that were not 
observed before. 



SUPERSYMMETRIC APPROACH 

Following the IWZ model^3], we consider the simplest 
case and divide our system into two boxes: one of the 
quantum dot (QD) and the other of the ring (see Fig.^l. 
Only the ring has a disorder given by a random GOE 
Hamiltonian -Hgoe of rank N and second momentum 
given by 
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FIG. 1: Aharonov-Bohm ring with an embedded quantum 
dot. 



2A is the radius of the Wigner semicircle. The QD is 
assumed to have one single level £d- We include the 
Aharonov-Bohm (AB) phase in the coupling term be- 
tween the walls of the QD and the ring. Doing so, the 
coupling is given by 
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Here vl,r are the coupling through the left and right 
side of the dot, which we can consider essentially the 
same. The AB phase is <fi = 7r$/$o, with $o = hc/e. In 
this way, the Hamiltonian of the system has dimension 
(1 + N) X (1 + N) and is given by 
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. V is a N dimensional vector with constant elements 

v{4>). 

The persistent current in the ring can be written in 
terms of the eigenenergies E n of the system as 
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with the condition that max£„ < Ep, Ep being the 
W Fermi energy. It is more convenient though to write this 
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expression using Green functions. In doing so we find 

(5) 

E' + and E'~ carry respectively a positive and a negative 
infinitesimal imaginary part and the limit is implicit. 

By performing the average over the impurities, we can 
not keep Ep fixed, since it changes both with impurity 



realization and the magnetic flux. The correct way to 
overcome this difficulty was proposed in Ref.jllj. and 
what we have to do is make an average over the num- 
ber of electrons too. This approach is ideal for treating 
an ensemble of rings, since due to uncertainties in their 
shapes, it is clear that the number of electrons is not the 
same in every ring, but must fluctuate about an average 
value N with an amplitude SN. This is equivalent to 
say that the Fermi energy is allowed to vary in an inter- 
val [Ei, Ei] of length S such that S = d5N, with d the 
mean level space at the Fermi surface. Thus we have the 
following expression for the current: 
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The bar denotes the average over the impurities. We 
see here that the current is expressed as a product of an 
advanced and a retarded Green function, which is essen- 
tially a two-point function. 

It is useful to define a generating functional Z from 
which we can obtain the current as simple derivatives, 
and furthermore, we can write our Hamiltonian in a 
Gaussian form such the procedure of averaging can be 
easily performed. We follow here the standard steps in 
the way presented in Ref. [J^- This generating function 
must then have the form 

Z[j,S(j)] = J d[^]cxp^^L^ 2 DL^ 2 ^ (7) 

where ^ is a graded vector of dimension 8(N + 1) and 
has the form: 
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and with the definitions 



D = (E--L-jk 2 )®l N+l 
V — v Q cos(0 + k\54>) (g) ei X 7v • 



e® 1 8 
V T H. 



V 

GOE C 



Here E = (E' + + E~)/2, u = E~ - E'+ and j and Scj) 
are source terms. L, ki and k 2 are 8x8 diagonal matrix 
and given by L = diag(l 4 , -1 4 ), ki = (1 2 , -I2, 0, 0), 
ki = (0,0,l2,—l2)- Finally l m represents the identity 
matrix of dimension mxm and ei x jv is a constant matrix 
of dimension 1 x JV and all elements equal 1. 

Hence the averaged persistent current is written in 
terms of the functional as 
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We may the proceed and perform the ensemble av- 
erage in a standard fashion followed by the Hubbard- 
Stratonovitch (HB) transformation. This yields for the 
generating functional 



Z = J d[cr]exp{£[cr]} 



where the Lagrangian is 
C[a] 

and 
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a is a 8 x 8 graded matrix which appears as an auxil- 
iary field to do the HB transformation. Trg denotes the 
graded trace, both over the 8x8 graded space and also 
the N + 1 degrees of freedom of the system. It is un- 
derstood here that the right bottom block of the above 
matrix has dimension 8iV x 87V. 

Since we are mostly interested in the limit N — > 00, 
we do the saddle point approximation to determine the 
main contribution for the exponent with u> = j = 5<j> — 0. 
It is useful to rewrite the logarithm of D as 
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where A = E - e— jL — jk 2 and T, N = (E 
jk 2 ) <8> ljv- Expanding the second term of Eq. (|13fl 
taking its trace, we can rewrite it as 



and 
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Now trg runs only over the graded space. 
We obtain for the saddle point equation 

A f 1 A% 2 cos 2 (0) \ 

° E-Xa\ N (E -e)(E- Act) - Nv 2 cos 2 (<f>) J ' 

(15) 

We note here that Nv 2 is of order of unity, and thus the 
second term in r.h.s. of Eq. I|15(l may be dropped in the 
limit N — ► oo. We find then the diagonal solution for the 
saddle point equation given by 

a u = ^-iAL, (16) 

where A = ^/l — (E/2X) 2 . This is not however the only 



solution. There is actually a manifold of solutions, which 
can be parametrized by some group generators T (see 
Ref. 0| to see its actual form) , and the general solution 
is 

<r S p = T-V D T = - iAQ . (17) 



We proceed with separation of a into the Goldstonc 
modes ctq an d the massive modes 8a, writing a = ctg + 
N~ 1 / 2 5cr and expanding the logarithm in function of 8a. 
We first note that uj, v 2 ~ N^ 1 , and then make the shift 
XaQ — > Actq — — jfc 2 , to find the effective Lagrangian 
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We have kept in the expansion only terms in first order 
in lu and j and up to second order in 8a. The linear terms 
in 8a coming from the logarithm cancel against those 
from a 2 . The integral over the massive modes is Gaussian 
and straightforward. Taking then the derivatives of the 
generating functional with respect to the sources j e 8<j) 
yields 



d 2 z 



djdScj) 



j=0=S^> 



trg[<r G fc 2 ] trg 



sA /sin(20) J d[<7 G ]e X p{£ } 

kiag + agki \ 

(£- e )-(7V Wo 2 cos 2 ^)/A) a G j • 



(19) 



where C Q = (Nlu/4X) trg[a G L}. 

Recalling Eq. I|17|l . the second trace of Eq. i|19fl can be 
put into a more convenient manner as 
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The change of variables that takes E and E' + into E 
and u implies a change in the region of integration, such 
that the new limits of integration are 
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Our model where we have used RMT to describe the 
disorder in the ring is valid for energies away from the 
boundaries of Wigner's semicircle. This implies that 
E <C A and therefore, we can disregard the term E/2X 
and take A = 1 in the integral over E. As a result, this 
integral becomes quite simple, and is equal to 
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Here, T((/>) = Nv 2 . cos 2 (0)/A. The averaged persistent 
current is then 
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SUPERSYMMETRIC INTEGRATION 

The integral f(uj) is done exactly with the parametriza- 
tion for Q given in Ref. and following the same steps. 
This is a lengthy and monotonous task, where we have 
to pay special attention to the anticomuttative variables. 
At the end of the day, we come to the final result for the 
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FIG. 2: Current as a function of the reduced flux "l>/$o, for 
coupling values To = Nvq/X. The current is given in units 
of cEf/^o- The parameters are E\ = 0,99 e E2 — 1,01, 
d — 10~ 3 , £d = 1,0, all in units of Ef- 



FIG. 3: Current as a function of the coupling To in units 
of cEf/^o- The parameters are E\ — 0, 99 e E2 — 1,01, 
d = 10 -3 , £d = 1,0, all in units of Ef- 



real part of f(u>) 



1G 



cos(27ru;/(i) — 1 
(ircu/d) 2 



- 16 



(25) 



We have used that the mean level spacing is given by 
d = nX/N. Defining T = Nv%/\, such that T{(/)) = 
Tq cos 2 (<f>) , we get the final expression for the current 
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RESULTS 

The plots of Fig. |2] to 0| show the current as function 
of different parameters as magnetic flux, resonance width 
To and energy of the QD level. The energy of the QD 
can be tuned experimentally through the gate potential 
V g . The natural energy unit adopted for the parameters 
is the Fermi energy of the ring Ep. For a semiconduc- 
tor ring constituted of a GaAlAs/GaAs heterostructure 
of mesoscopic dimensions, i.e., with a diameter of fj,m, 
the total electron number is of thousands and the Fermi 
energy laying in the meV. The electron effective mass 
is about 0.07 the free electron mass[l^| and the elastic 
mean free path is about the circumference of the ring. 



In the first plot we observe that the periodicity of the 
current is <&o, double of the value found in the case with- 
out the QD. First, we should remember that the per- 
sistent current in a single isolated ring has a periodicity 
with the magnetic flux of $0 both for the case of an odd 
or even number of electrons. What differs the current in 
each case is just a shift of $o/2- Thus, due to the very 
symmetry of those plots, when we perform an average 
over the number of electrons, i.e., we sum the currents 
in both cases, we end up with a current with a period- 
icity of <£n/2. This is valid for a ring with or without 
disorder [l5|. The conclusion we can derive here is that 
the very presence of the QD breaks down the symmetry 
of the current for an even or odd number of electrons. 
Such asymmetry was also observed for a clean ring in the 
Kondo regime|j|, but now we can see that this feature is 
even more general. 

Another remark to be done is about the coupling con- 
stant. For small values of this quantity (Tq < 10~ 3 Ep), 
the current is basically a sine and it scales with Tq, such 
that I /Tq is an universal function. For growing values of 
Tq the current form begins to alter, the maximum (and 
minimum) of the curve approaches $o/2 and its ampli- 
tudes grows a bit slower with vT . In Fig. E|we see how 
the current behaves as a function of Tq for some selected 
values of the magnetic flux. The first and obvious con- 
clusion is that the current vanishes when Tq — > 0, then 
this case is like we had opened the ring and no current 
is allowed to flow. In the same manner for Tq very big, 
what we get again is a vanishing current. What happens 
now is something different. The life time of the reso- 
nance is inversely proportional to Tq and therefore the 
probability of the electron to pass through the QD is di- 
minished. If there is experimentally a controllable way 
to alter the coupling between the ring and the dot, then 
this plot shows where we can get a maximum amplitude 
for the current. 

Our main result, however, may be observed in Fig. 0] 
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FIG. 4: Current as a function of QD energy ed in units of 
cEf/$o- The parameters are E\ = 0, 99 e E2 — 1,01, d — 
10" 3 , To = 0, 1, all in units of E F , and $ = 0, 3$ - 



where we have analyzed the variation of the current am- 
plitude as function of the energy of QD level. We have 
observed a peak in the current for ed near the Fermi en- 
ergy of the ring. Such result is expected if we note that, 
in the case without the QD, the main contribution for the 
total current comes from the current associated with the 
electron in the uppermost level. Thus, when the energy 
of the QD coincides with this energy, there is a peak in 
the current indicating the resonance between both levels. 
The second peak however is quite curious. It occurs at 
about half of the Fermi energy and is symmetric to the 
first, only with opposite signal. The explanation for this 
phenomenon is the following. The coupling between the 
energy levels in the ring and that of the QD gives then 
a width and a Lorentzian shape. Since the current is es- 
sentially a derivative of the energy with respect to the 
magnetic flux, we get that the individual level contribu- 
tion for the total current is proportional to the difference 
of its energy and ed- Thus, when the QD level crosses 
down the Fermi energy, the contribution of the levels 
with higher energies than ed have opposite signal. For 
a given value of ed, the total current vanish. Lowering ed 
even more, the current changes signal and increases until 
it reaches a maximum value and then decreases again as 
ed is much lower than the energy levels of the ring. 

Unlikely the results obtained in other works[2, La-EIEI- 
this extra peak is neither related to the Fano effect nor 
to the Kondo effect. The Fano effect HI El El occurs 
when there is an interference between two paths through 
which the electron can pass, one with discrete levels, and 
the other a continuous band. To observe such effect, one 



may use one ring with a QD and two leads connected in 
opposite sides of the ring[g. Measuring the conductance 
we find peaks with asymmetric line shapes. On the other 
hand, Kondo effect 0,0,0] is related to the fact that the 
QD behaves like a magnetic impurity in the ring, and to 
observe it, it is necessary that the net spin in the dot be 
different of zero. As a result of the correlations between 
the electrons in the QD and those of the conduction band, 
there is a resonance near the Fermi energy with works as 
an extra channel for the electron tunneling. Our model 
does not account for the electronic spin and therefore 
there is no way of interpreting our result as a Kondo 
resonance, furthermore the second peak lays far below 
the Fermi level. In conclusion, the second peak arises 
as a constructive interference phenomenon between the 
levels of the ring and is purely a manifestation of the 
quantum nature of the system. 
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